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1.  Introduction 


This  report  summarizes  the  scientific  and  technical  progress  made  during  the  12-month 
duration  May  1,  1990  through  May  31,  1991  under  the  subject  contract. 

Since  this  is  the  initieJ  year  of  the  contract  and  CRRES  science  data  usable  for  mod¬ 
eling  were  not  available  until  about  September  1990.  our  modeling  effort  up  to  the  present 
is  limited  to  infrzistructure  and  single  cases  of  dynamic  storm  effects,  as  specified  in  the 
original  proposal.  The  major  efforts  on  a  storm  time  model  with  statistical  support  await 
the  formation  of  an  outer  belt  storm  data  base,  which  is  being  accumulated  currently.  In 
spite  of  the  short  time  interval  since  initiation,  at  least  four  significantly  new  results  in  ba¬ 
sic  magnetospheric  science  have  been  obtained  under  the  auspices  of  this  contract.  These 
are:  (1)  a  complete  and  exact  solution  of  the  outer  magnetosphere  with  properties  that 
mimic  the  bow  shock,  the  magnetosheath,  the  magnetopause  and  IMF  merging  for  any 
given  shape  of  the  magnetopause  and  botv  shock  boundaries;  (2)  discovery  of  the  general 
properties  of  high-energy  (-^1  MeV)  electron  and  proton  encounters  with  the  outer  magne¬ 
tospheric  structures,  leading  to  a  comuputer  verification  of  the  concept  of  outer  radiation 
belt  particle  recycling  by  geomagnetic  storms:  (3)  new  interpretations  for  the  causes  of  the 
violent  fluctuations  of  outer  belt  electron  fluxes  as  results  of  insights  gained  in  (1)  and  (2) 
above;  and  (4)  a  new  integral  view  of  the  initial  and  recovery  phaises  of  outer  belt  dynamic 
responses  to  geomagnetic  storms.  The  leist  item  is  the  result  of  insight  from  our  current 
examination  of  CRRES  data  in  the  context  of  our  theoretical  development.  Although  it 
is  too  early  to  determine  if  the  modeling  strategies  developed  hither  to  may  or  may  not 
require  major  modification,  as  result  of  compiling  a  sufficient  stormtime  data  base;  it  is 
significant  to  note  that,  up  until  this  effort  organized  by  the  CRRES  program,  there  has 
not  been  an  attempt  to  develop  an  integrated  model  of  the  dynamiced  response  of  the  outer 
electron  belt  in  the  context  of  magnetic  storms,  in  a  quantitative  manner  according  to  the 
response  times  of  its  phases. 

In  this  report,  we  shall  attempt  to  give  a  scientific  overview  of  the  current  modeling 
strategy,  rationale  and  observational  motivation,  as  if  the  theoretical  approaches  are  well 
established.  This  by  no  means  implies  that  we  shall  not  modify  the  modelling  approaches 
in  response  to  the  accumulating  CRRES  data  base.  Rather,  we  believe  that,  unless  we 
commit  ourselves  to  the  test  of  a  specific  organized  conceptual  modelling  scheme,  the 
well-established  violent  fluctuations  of  the  outer  belt  will  not  be  easily  unravelled. 

The  main  body  of  this  report  is  divided  into  6  sections:  (2)  Outer  Magnetosphere 


Modeling,  for  establishing  an  infrastructured  context  for  viewing  the  outer  belts;  (3)  Dif¬ 
fusion  Modeling,  for  studying  responses  on  the  time  scale  of  days;  (4)  Particle  Encounters 
with  Outer  Magnetosphere  Structures,  for  studying  responses  on  the  drift  time  scale;  (5) 
Test  Case;  (6)  Future  Work;  and  (7)  Publication  and  Presentation  List. 

2.  Outer  Magnetosphere  Modeling 

To  determine  the  response  of  the  outer  belt,  it  is  necessary  to  have  a  model  of  the 
outer  magnetosphere  that  mimics  the  magnetospheric  response  to  solar  wind  parameters. 
An  outer  magnetosphere  model  must  include  the  responses  of  the  outer  magnetosphere 
structures  such  as  bowsliock,  magnetopause,  magnetosheath  and  merging  efficiency  pa- 
rameterizations.  There  is  currently  no  such  basic  physics  model  as  the  detailed  plasma 
physics  axe  still  being  determined.  There  are  magnetohydrodynamic  and  potential  field  so¬ 
lutions  which  represent  the  current  state-of-the-art;  all  of  these  models  require  CRAY-like 
computation  requirements  with  long  running  times.  Such  models  are  scientifically  signifi¬ 
cant  but  entirely  impractical  in  the  context  of  usage  for  outer  belt  dynamic  modeling  bcised 
on  CRRES  data;  since  the  role  of  the  outer  magnetosphere  model  is  to  serve  ais  part  of 
the  infrastruction  for  the  dynamical  motion  of  energetic  particles.  On  the  other  hcind,  the 
reality  is  that,  without  a  reasonable  analog  of  ihe  outer  magnetosphere  structure,  espe¬ 
cially  on  the  dayside,  little  organized  progress  can  be  made  in  the  context  of  CRRES  data 
interpretation. 

Stem  (JGR,  1986)  described  an  elegant  approach  to  the  magnetosphere  modeling 
problem  by  demonstrating  that  the  difficult  potential  boundary  value  problem  can  be  sim¬ 
plified  by  choosing  a  co-ordinate  system  in  which  the  magnetospheric  boundarj-  structures, 
such  eis  bow  shocks  or  magnetopauses,  lie  on  a  co-ordinate  surface.  He  demonstrated  how 
an  elegant  inner  magnetosphere  solution  with  appropriate  currents  to  yield  a  paraboloidal 
magnetopause  can  be  specified  efficiently  with  an  expansion  of  the  paraboloidal  eigenfunc¬ 
tions  for  Laplace's  potential  equation. 

In  the  first  phase  of  our  outer  magnetosphere  modeling  effort,  we  adopted  this  ap¬ 
proach  to  obtain  exact  solutions  that  include  paraboloidal  outer  magnetospheric  structures 
(bow  shock  and  magnetopause),  including  dayside  merging  of  the  IMF  and  the  internal 
geomagnetic  field  represented  by  an  efficiency  parameter.  The  price  paid  for  such  an  el¬ 
egant  solution  is  that  the  shapes  of  outer  magnetospheric  structures  are  restricted  to  be 
on  a  co-ordinate  surface;  and,  in  addition,  tail  effects  are  ignored.  The  latter  drawback 
is  not  serious  for  outer  belt  modeling  because  the  outer  belt  does  not  extend  to  £-values 
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much  higher  than  8.  Figure  1  summarizes  the  mathei.',aticE  involved.  Figure  2  shows  a 
typical  sample  of  the  exact  potential  solution  result  that  is  virtually  indistinguishable  from 
magneiohydrodynamic  calculations  requiring  extensive  computational  resources.  We  have 
used  this  paraboloidal  ouier  magnetosphere  code  to  mimic  dayside  outer  magnetospheric 
structures  for  various  test-Ccise  scenarios  which  we  will  discuss  in  Section  4  and  5. 

The  paraboloidal  outer  magnetosphere  model  has  the  virtues  of  simplicity  and  time¬ 
liness  for  study  of  the  first  results  of  CRRES  outer  radiation  belt  data.  Its  elegance  as 
an  exact  solution  hzLS  intrinsic  scientific  vEdue;  however,  two  major  elements  need  to  be 
accommodated;  (1)  a  magnetosphere  inside  the  magnetopause  with  a  tail  structure  and 
a  non-paraboloidal  magnetopause  shape  w'hich  is  consistent  with  observations  and  with 
the  existence  of  a  tail.  While  there  are  several  such  magnetosphere  models  available  to 
be  accommodated  into  our  development  of  the  outer  magnetosphere  boundary  conditions, 
the  development  of  the  Source  Surface  Model  (SSM)  prescription  by  Schulz  and  McKab 
(GRL  1989)  has  the  virtue  that  it  can  accommodate  arbitrary  magnetopause  shapes  ai.d 
their  appropriate  tail  currents.  We  have  recently  initiated  a  merger  of  the  Outer  Magne¬ 
tosphere  Potential  Model  with  the  Source  Surface  Model.  The  results  of  this  merger  are 
very  satisfactory. 

Figures  3  and  4  show  such  implementations  for  a  paraboloidal  bow  shock  and  a  magne¬ 
topause  shape  that  fit  observations  well.  Only  outer  magnetosphere  field  lines  are  shown  in 
these  cases  for  simplicity.  Computation  requirements  for  this  type  of  model  are  about  one 
order  of  magnitude  higher  than  the  paraboloidal  model,  however,  it  is  still  economical  com¬ 
pared  to  magnetohydrodynamic  solutions.  This  work  is  in  collaboration  with  Mike  Schulz, 
who  is  not  supported  by  the  subject  contract.  An  invited  paper  on  this  significantly  new 
development  in  basic  magnetosphere  physics  will  be  presented  by  Mike  Schulz  in  the  lAGA 
meeting  in  Vienna  in  .4ugust  1991.  A  detailed  paper,  under  the  leadership  of  Lockheed 
and  the  sponsorship  of  this  subject  contract,  is  being  prepared  for  publication.  Because 
of  the  heavy  computation  requirements,  development  of  this  model  is  partially  funded  by 
the  Lockheed  Independent  Research  Program  in  addition  to  the  subject  contract. 

3.  Diffusion  Model’ng 

A  reasonable  outer  magnetosphere  model  that  mimics  its  response  to  solar  wind  con¬ 
ditions  is  an  essential  infrastructure  element  in  dynamic  modeling  of  the  outer  bell.  The 
other  essential  element  is  an  efficient  accounting  of  the  particle  diffusion  processes  that  are 
expected  to  take  place  in  the  re-establishment  of  the  equilibrium  after  the  initial  storm 
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response. 

In  the  fore-runner  study  with  SCATHA/SEP  data,  we  demonstrated  that  simultane¬ 
ous  radial  and  pitch  angle  diffusions  can  accoimt  for  the  behavior  of  outer  electron  belt 
dynamics  on  the  time  scale  of  a  day  or  so.  The  physical  diffusion  parameters  derived  from 
the  SCATHA  high  energy  electron  data  are  in  good  agreement  with  other  independent 
determinations  and  theoretical  expectations.  For  CRRES  dynamic  modeling,  we  make  use 
of  this  facility  that  has  been  constructed  in  the  last  three  years.  The  principles  of  the 
diffusion  modeling  have  been  published  in  two  detail  papers  [JGR  93,  2619,  1988,  JGR  95, 
12069, 1990]. 

/.Ithough  the  basic  assumptions  and  applications  of  the  simultaneous  radial  and  pitch 
angle  diffusion  model  are  the  same  between  the  SCATHA  tests  and  the  CRRES  applica¬ 
tions,  there  is  still  a  major  difference  that  needs  to  be  addressed.  The  cases  treated  in  the 
SCATHA  data  base  are  in  moderate  geomagnetic  activity  without  consideration  of  drift- 
time  scale  responses  to  clearly  defined  geomagnetic  storms;  whereas,  the  main  objective 
of  the  CRRES  storm  response  model  is  to  attack  the  outer  belt  response  to  storms  at  ail 
time  scales. 

Aside  from  issues  such  as  the  geomagnetic  f.eld,  the  diffusion  response  in  the  recovery’ 
phase  of  the  response  model  requires  an  additional  term  on  the  right  hand  side  of  the 
diffusion  equation;  namely, 

a/  ... 

dt~  '  '  r 

where  A  is  the  simultaneous  radial  and  pitch  angle  diffusion  operator  specifying  the  time 
rate  of  phase  space  density  change  (-1^)  due  to  diffusive  processes.  The  new  term,  with 
time  seeding  7,  represents  a  ’  source"  term  due  to  non-diffusive  time  rate  of  change  of  phase 
space  density  {F).  In  this  regard,  the  nomenclature  ’’source’’  is  a  misnomer  originating 
from  tradition,  which  sweeps  anything  unaccounted  for  by  diffusion  as  a  ’’source’*.  In  our 
storm  time  response  model,  for  example,  parts  of  this  soiurce  term  can  be  the  population 
of  outer  belt  electrons  removed  to  the  outer  mangetosphere  by  non-diffusive  processes 
such  as  scatterings  on  drastically  changed  magnetospheric  structures.  The  recycling  of 
such  populations  into  the  usual  outer  belt  location  by  radial  diffusion  during  the  recovery 
response  of  the  outer  belt  is  represented  in  the  diffusion  equation  as  a  particle  source  of 
non-diffusion  origin.  The  objectives  of  a  dynamic  modeling  effort,  ciside  from  deriving 
diffusion  parameters  (properties  of  A),  is  to  specify  the  ’’source"  term  that  would  account 
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fur  the  properties  of  a  CURES  stoi  i  data  bsise,  which  is  currently  being  accumulated. 
At  the  present  time,  several  interesting  aspects  of  F  have  been  revealed  by  examining  the 
detailed  phase  space  density  of  a  test  case  (storm  of  August  1990),  which  is  discussed  in 
Section  5  in  more  detail. 

Meanwhile,  the  computational  procedure  to  unravel  the  diffusion  equation  (1)  to  spec¬ 
ify  the  diffusion  parameters  of  A  and  to  construct  the  appropriate  F  requires  decomposition 
of  the  observed  distribution  function  /  into  eigenfunctions  of  A  so  as  to  facilitate  extrac¬ 
tion  of  the  diffusion  constants.  Since  we  have  chosen  an  approximate  representation  of  A 
such  that  its  eigenfunctions  are  also  a  complete  set  in  pitch-angle  space,  the  rendering  of 
CRRES  data  in  terms  of  (1)  requires  the  first  step  of  decomposing  /  into  pitch-angle  eigen¬ 
functions  that  allows  us  to  reduce  the  independent  dimensions  of  (1)  by  one.  The  critical 
test  of  this  step  is  whether  the  pitch  angle  decomposition  of  /  into  approximate  pitch-angle 
eigenfunctions  (Bessel  functions  of  order  zero)  is  as  efficient  as  the  SCATHA  test  cases. 
It  is  to  noted  that  the  CRRES  modeling  case  requires  this  decomposition  throughout  the 
storm  phases.  We  demonstrated  this  decomposition  in  Figures  5  and  6. 

In  Figure  5,  the  pitch-angle  decomposition  amphtudes  Qn  in 


/  =  E<2" 


^nO  (1  ~  ^no) 


y/2Jo{k„xlXf.) 


(2) 


is  shown  for  the  entire  CRRES  orbits  of  77  (initial  storm  response)  and  82  (late  recovery 
response)  for  electron  energy  channels  800  keV  and  400  keV  respectively.  The  dependence 
of  these  amplitudes  Qn  upon  L  need  not  concern  us  for  the  moment.  The  important  point 
is  that  in  the  outer  belt  the  amplitudes  of  the  n  =  0  -  4  modes  ccount  for  the  total  of  (2), 
even  for  the  violently  changing  initial  response  (orbit  77). 

We  have  completed  the  pitch  angle  decomposition  of  the  data  base  (orbits  76-82)  for 
the  sample  storm  of  August  1990  and  foimd  the  decomposition  to  be  rapidly  and  stably 
decreasing  with  order.  .4s  in  the  SCATHA  tests,  4  orders  of  pitch  angle  eigenfunctions  are 
quite  sufficient  to  decompose  the  CRRES  test  case  above. 


4.  Particle  Encounters  with  Outer  Magnetosphere  Structures 

.4  third  infrastructural  element  is  required  to  examine  the  behavior  of  the  outer  belt 
in  the  initial  response  to  geomagnetic  storm  conditions.  The  fact  that  the  dayside  magne¬ 
tosphere  is  highly  distorted  in  the  storm  initial  phaise  is  evident  and  discussed  in  Section 
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2.  Of  chief  interest  to  outer  belt  dynamic  modeling,  however,  is  how  do  the  trapped  en* 
ergetic  particles  move  in  the  distorted  outer  magnetosphere  on  the  drift  and  bounce  time 
scales?  The  answer  to  this  question  is  crucial  to  the  description  of  the  initieil  pheise  of  a 
dynamic  response  model  and  to  construction  of  the  non-diffusive  "source”  function  F  for 
the  recover}'  phase. 

Crucicil  cis  the  description  of  energetic  particle  encounters  with  the  outer  magneto- 
spheric  structures  may  be,  it  is  also  unfortunate  that  there  is  no  economical  method  to 
study  the  properties  of  /,  as  an  ensemble,  except  to  exzunine  and  generalize  the  behavior 
of  individual  particle  traces  in  terms  of  a  simulation.  We  have  done  extensive  studies  in 
this  area  at  the  early  part  of  this  contract  year  while  awaiting  CRRES  data  to  accumulate. 

We  have  accumulated  an  extensive  computation  data  base  of  MeV  electron  and  pro¬ 
ton  trajectories  of  various  initial  pitch  angles  moving  in  paraboloidal  outer  magnetospheres 
such  as  those  shown  in  the  next  section.  In  addition,  we  have  also  included  outer  magne¬ 
tospheres  with  superimposed  large  scale  wave  structures  in  the  computational  trajectory- 
data  base.  This  large  data  base  of  trajectories  cannot  be  conveniently  exhibited  here  be¬ 
cause  of  variations  in  initial  pitch  angle  and  outer  magnetosphere  geometiy;  therefore,  the 
samples  shown  in  the  next  section  for  specific  outer  magnetospheres  are  more  appropriate 
samples  to  indicate  the  inventory  of  our  study  results  in  this  area.  Here,  we  state  several 
general  conclusions  that  are  applicable  to  MeV  particle  motion  even  outside  of  the  context 
of  outer  belt  modeling: 

(1)  MeV  particle  encounters  with  outer  magnetosphere  structure,  as  simulated  fav¬ 
our  parabolodial  magnetopause,  magnetosheath  and  bow  shock,  rarely  lead  to  loss  of  the 
particle  into  the  dayside  heliosphere. 

(2)  The  most  common  result  of  such  encounters  leads  to  scatter  of  the  particle  trajec¬ 
tory  from  tightly  defined  trapped  bounce-drift  motion  to  quasi-bounce-drift  motion  steered 
by  the  outer  magnetosphere  structure  on  trajectories  that  span  a  much  Icirger  volume  of 
space  but  mostly  inside  the  magnetosphere.  The  immediate  consequence  of  this  feature  is 
the  dispersal  of  outer  belt  fluxes  to  a  low  level  in  those  populations  of  the  outer  belt  that 
have  encountered  the  outer  magnetosphere. 

(3)  Electrons  (protons)  invariably  enter  the  magnetosphere  on  the  dawn  side  (dusk 
side)  and  leave  on  the  dusk  side  (dawn  side). 

(4)  Since  outer  magnetosphere  models  used  in  this  series  of  studies  do  not  contain  a 
tail,  the  eitects  of  the  geomagnetic  tail  structure  to  MeV  particles  is  not  revealed  here. 
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Other  studies  at  lower  energy'  also  indicate  scattering  but.  again,  contained  inside  the 
magnetosphere. 

The  reader  is  referred  to  Section  5  for  detailed  ramification  of  these  conclusions  in 
terms  of  the  specific  test  case  of  the  August  1990  storm. 

5.  Test  Case:  Storm  of  August  1990 

Despite  the  small  magnitude  of  geomagnetic  actinty  associated  with  the  storm  of  Au¬ 
gust  1990,  its  effects  on  the  outer  magnetosphere  and  the  advantageous  location  of  CRRES 
on  the  momside  side  flanks  of  the  magnetosphere  combined  to  yield  distinctive  features 
that,  in  our  opinion,  have  sharpened  our  understanding  of  the  outer  belt  dynamic  response 
to  solar  wind  conditions.  This  section,  as  well  as  the  major  pari  of  our  interpretation  and 
modeling  effort  in  the  latter  half  of  the  contract  period,  is  devoted  to  this  storm,  not  only 
£LS  a  sample  of  our  standard  modeling  procedure  that  will  be  applied  to  every  storm  of  our 
accumulating  CRRES  Storm  Data  Base,  but  also  as  a  stand-alone  science  ctnadysis  project 
since  we  are  able  to  isolate  several  new  features  of  the  outer  belt  response.  The  first  part 
of  this  section  deals  with  the  organized  SEP  data  bsise  for  this  event  while  the  second  pzul 
deals  with  modeling  the  initial  phase  of  the  outer  belt  dynamic  response.  The  recovery 
phase  of  the  response  has  yet  to  be  modeled. 

a.  SEP  Data  for  the  Event 

The  solar  wind  conditions  for  the  event  is  capsulized  in  Figure  7  which  shows  the 
subsolar  standoff  distances  of  the  magnetopause  and  the  IMF  components  as  the  initial 
phase  of  the  event  progresses.  The  standoff  distance  shows  two  sudden  inward  motions.  We 
have  used  the  information  in  this  figure  to  define  analogous  outer  magnetosphere  models 
for  both  inward  motions  of  the  standof  andf  distance.  Only  results  of  the  second  inward 
motion  (to  L  =  6.5)  will  be  discussed  here  since  immediate  responses  of  the  outer  belt 
energetic  electrons  are  observed  by  CRRES  in  its  apogee  region  shown  on  Figure  8.  The 
electron  flux  spectographs.  Figure  9,  show  the  electron  fluxes  of  the  SEP  instruments  for 
orbits  75  (before  the  event),  77  and  79,  ordered  by  L-value.  Focusing  upon  the  energetic 
band  (above  200  keV).  it  is  quite  clear  that  a  "drop-out"  of  flux  was  detected  within  the 
drift  time  scale  (less  or  order  of  one  hour).  Recovery  to  pre-event  fluxes  for  these  particles 
had  not  occurred  one  day  later.  Indeed,  it  took  2-3  days  before  complete  flux  recovery. 
The  short  time  scale  of  the  flirx  drop  out  is  the  first  notable  feature  because  it  clearly  do^ 
not  indicate  a  diffusion  process. 


Next,  the  pitch  angle  characteristics  of  the  energetic  electrons  cilso  yield  significant 
insight  into  the  dynamics  of  the  dropout.  Figures  10-13  is  a  set  of  pitch  angle  distributions 
measured  on  orbit  77  at  two  L- values  before  (inbound)  and  after  (outbound)  the  flux  drop 
out  mapped  to  the  equator  for  uniformity  of  comparison,  as  indicated  by  the  flux  levels  for 
both  the  400  keV  and  800  keV  channels.  Two  important  features  are  apparent  when  the 
set  is  examined  as  a  whole.  First,  the  inbound  (before)  pitch-angle  distribution  exhibits 
the  usual  pitch-angle  diffusion  signature  sloping  towards  smaller  pitch  angles;  whereas  all 
the  outbound  (after)  distributions  are  isotropic  up  to  the  loss  cone.  The  dynamical  process 
responsible  for  the  flux  drop  out  must  be  able  to  isotropize  the  pitch-angle  distribution 
within  a  time  less  than  3  hours  shown  on  the  L  =  6.5  pair,  Figures  12  and  13.  That  the  flux 
drop  out  process  is  temporal  is  also  clearly  indicated  by  the  examining  the  magnetic  local 
time  (MIT)  differences  of  the  inbound  and  outbound  fluxes.  The  ’’before”  and  ’’after” 
states  art  clearly  marked  within  a  distance  of  1.5  hour  MLT.  A  diffusive  or  convective 
process  that  can  produce  a  Aua  decrease  of  800  keV  electrons  of  up  to  an  order  of  magnitude 
within  a  spatial  extend  of  1.5  hour  MLT  is  difficult  to  imagine.  We  suggest  a  scattering 
process  such  as  described  in  the  previous  section  and  analogously  modeled  in  the  latter 
part  of  this  section. 

Another  important  characteristic  of  the  outer  belt  electrons  above  400  keV  before, 
during  and  after  the  flux  drop  out  is  revealed  in  the  set  of  Figures  14,  15  and  16.  These 
spectra  in  the  energy  interval  400  keV  to  about  1.5  MeV  exhibit  the  same  decrease  of 
roughly  1.5  decade  of  flux.  Apparently,  the  dynamics  of  the  flux  drop  out  does  not  effect 
the  spectrum  in  this  energy  range  which  is  the  concern  of  our  model.  The  level  of  multi- 
MeV  electrons  also  seem  to  remain  constant,  i.e.  does  not  participate  in  the  drop  out, 
and  the  spectra  of  the  lev/  energy  electrons  art  drastically  different  in  these  Figures.  The 
latter  can  be  due  to  dB/dt  (inductive  electric  field  effects),  but  the  unchanged  spectra  of 
the  high  energy  electrons  do  not  fit  the  inductive  electric  acceleration  hypothesis  as  the 
amount  of  acceleration  is  proportional  to  the  initial  energy  of  the  particle,  thus  resulting 
in  spectrum  changes. 

The  only  hypothesis  that  agrees  with  all  of  the  above  observed  features  that  we  have 
been  able  to  come  up  with  is  that  of  scattering  of  the  outer  belt  population  giving  rise 
to  an  increased  volume  in  which  the  population  is  contained.  The  scatterings  with  outer 
magnetosphere  features  can  reasonably  isotropize  the  pitch  angle  distribution  in  a  short 
time  without  affecting  the  spectrum. 
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b.  Analogy  Modeling  for  the  Event 

In  this  segment,  the  entire  effort  was  focused  on  analogy  modeling  for  the  initial  part 
of  the  outer  electron  belt  and  trapped  proton  response  to  the  storm  of  CRRES  orbit  77. 
During  this  period,  the  outer  magnetosphere  is  very  dynamic.  The  dynamic  response  time 
scale  of  particle  distributions  is  considerably  shorter  than  normal  diffusion  times;  therefore, 
a  diffusive  picture  is  not  appropriate  for  these  initial  orbits.  This  view  is  clearly  supported 
by  observations  cited  in  the  previous  segment. 

The  focus  of  the  modeling  effort  is  to  present  a  clear  picture,  for  the  first  time  as  we 
can  determine  to  date,  of  the  initieJ  response  of  the  outer  radiation  belt  to  sudden  changes 
of  solar  wind  flow  and  magnetic  parameters.  We  are  using  the  storm  as  a  benchmark  for 
our  outer-belt  dynamic  modeling  work.  At  the  present  time,  sufficient  details  of  our  data 
analysis  and  parsJlel  simulation  work  have  emerged  to  afford  a  strikingly  exciting  and  ne;v 
view  of  radiation  belt  processes  occurring  in  the  outer  magnetosphere.  We  shall  present 
a  series  of  charts  here  to  document  the  current  state  of  several  parallel  segments  of  the 
initial  storm  response  modeling  effort.  The  segments  of  data  and  modeling  are  shown  in  a 
logical  sequence,  as  follows: 

1.  Solar  wind  stand-off  distance  and  IMF  data  2.  Outer  magnetosphere  model  with 
paraboloid  magnetopause  and  bow-shock  corresponding  to  the  solar  wind  data.  3.  The 
storm  response  data  indicating  an  extensive  ’’drop  out”  of  outer  belt  flux.  4.  Study  of 
radiation  belt  electron  and  proton  motion  in  the  region  affected  by  the  sudden  movement 
of  the  magnetopause  inward. 

Taking  averages  of  the  standoff  distance  and  IMF  components  for  a  period  of  about 
3  hours  before  and  after  the  second  sudden  change  of  Figure  7,  we  arrived  at  critical 
parameters  to  model  the  outer  regions  of  the  magnetosphere.  Figure  17  shows  the  X-Z 
plane  projection  of  our  open  outer  magnetosphere  model  with  merging  efficiency  of  0.2 
from  the  period  before  the  sudden  compression.  The  IMF  component  values  are  shown  at 
the  top.  The  outer  magnetosphere  model  is  an  exact  solution  of  Laplace’s  equation  with  a 
paraboloidal  shape  that  fits  the  observed  average  stand  off  distance  and  IMF  magnitude 
parameters. 

Figures  18, 19  and  20  give  the  three  orthogonal  projections  of  the  outer  magnet ospheric 
field  after  the  sudden  compression.  The  IMF  components  and  other  parameters  are  given 
at  the  top.  Comparisons  of  Figures  17  and  18,  19  and  20  show  the  drastic  changes  of  the 
outer  magnetosphere  configurations. 
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To  model  what  happens  to  radiation  belt  particles  trapped  in  the  region  of  dynamic 
change,  the  vicinity  of  6-8  RE  as  measured  on  the  dayside,  we  investigate  the  dynamical 
behavior  of  such  particles  (10  MeV  electrons  and  protons  for  the  sake  of  reducing  computer 
running  time)  before,  during  and  after  the  outer  magnetosphere  configuration  change. 

We  have  made  many  runs  of  electron  and  proton  trajectories.  The  following  is  just  a 
sample  that  brackets  0  and  90  degree  equatorial  pitch  angles.  The  trajectories  were  started 
at  90  degree  pitch  angle  near  the  mirror  point  of  field  lines  that  were  just  inside  of  the 
magnetopause  to  assure  ’’classical”  trapping  or  ’’scattered”  trapping  (’’rattling  along  the 
magnetopause  inside”)  that  we  have  discussed  in  previous  reports.  Depending  on  minor 
differences  of  the  starting  pitch  angle,  the  particles  ”fen  out”  into  equatorial  pitch  angles 
of  0-90  degrees.  The  parameters  of  each  set  of  three  projections  per  sample  trajectory  are 
printed  at  the  top  of  each  figure. 

Figures  21,  22  cind  23  show  the  projections  of  an  electron  trapped  in  the  ’’Before” 
outer  magnetosphere  with  equatorial  pitch  angle  of  about  10  degrees.  This  particle  is  very 
near  the  ’’Before”  loss  cone.  Figures  24,  25,  26  show  the  same  projections  of  the  same 
particle  in  the  ’’After”  outer  magnetosphere.  It  is  seen  that  it  did  not  bounce  even  once 
before  it  went  off  into  the  tail  region.  Since  our  outer  magnetosphere  model  is  not  expected 
to  be  good  in  the  tail,  we  cannot  comment  on  the  fate  of  this  particle. 

Figures  27,  28,  29  depict  projections  of  an  electron  with  equatorial  pitch  angle  near  90 
degrees  in  the  vicinity  of  6RE  near  the  morning  dayside  in  the  ’’Before”  magnetosphere. 
Figures  30,  31,  32  depicts  the  same  particle  in  the  ’’After”  magnetosphere.  It  is  seen 
that  this  particle,  in  the  vicinity  of  CRRES,  is  scattered  into  the  quasi-trapped  trajectory 
’’rattling”  near  the  inside  wall  of  the  magnetopause.  Considering  that  this  particle  has  90 
degree  equatorial  pitch  angle  before,  the  ’’drop  out”  of  flux  in  the  L=6  region  shown  in 
Figure  9  is  easily  understood.  The  important  point,  however,  is  that  this  particle  is  not 
lost  to  the  magnetosphere.  It  is  merely  dispersed  into  the  ’’rattling”  trajectory.  Later, 
perhaps  much  later  as  our  flux  drop  out  did  not  get  filled  in  completely  even  after  Orbit 
84,  this  particle  can  be  recycled  back  into  the  outer  belt  by  diffusive  processes. 

The  next  set  of  figures  investigate  the  effects  of  the  sudden  decrease  of  the  magne¬ 
topause  standoff  distance  upon  a  stably  drifting  particle.  Because  of  the  long  computation 
time  for  electrons,  we  use  a  10  MeV  proton  as  test  particle.  Figures  33,  34,  and  35  shows 
the  stably  drifting  trajectory  of  the  proton  in  the  ’’Before”  outer  magnetosphere.  Note 
that  this  particle  has  equatorial  pitch  angle  in  the  45  degree  region.  After  the  decrease  of 
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the  standoff  distance,  shown  in  Figures  36,  37,  and  38  we  see  that  the  stable  drift  persists 
except  that  it  now  ’’rattles”  along  the  magnetopause  on  the  dayside,  staying  on  some  of 
the  time  in  the  magnetosheath. 

Conclusion:  We  have  made  i.  great  deal  of  progress  understanding  the  initial  response 
of  the  outer  belt  to  magnetic  storms  by  modeling  both  the  magnetic  field  and  particle 
behavior  of  the  selected  CRRES  event.  Solar  wind  data  is  absolutely  necessary  to  achieve 
this.  This  flux  ’’drop  out”  at  the  initial  phase  of  the  magnetic  storms  is  likely  to  be 
understood  as  a  dispersal  of  the  trapped  outer  belt  by  scattering  with  outer  magnetospheric 
structures.  Diffu.3ive  processes  operate  over  much  longer  time  scalds  later  to  fill  in  the  void 
left  by  scattering. 

c.  Recovery  Phase 

As  indicated  previously,  the  diffusion  model  for  the  recovery  phase  of  the  event  is 
currently  underwa}'.  Figures  39  and  40  show  the  pitch-angle  distributions  typical  of  diffu¬ 
sion.  Further,  these  figures  with  fluxes  from  inbound,  outbound  and  the  two  halves  of  the 
magnetic  field  lines  mapped  to  the  equator  show  imequivocally  that  there  are  no  local-time 
or  temporally  sudden  processes  in  operation. 

6.  Future  Work 

In  the  above,  we  have  outlined  the  salient  observational  features  of  the  CRRES /SEP 
outer  belt  dynamic  data  especially  in  the  context  of  storm  responses.  Further,  we  have 
put  together  a  plausible  modeling  hypothesis  to  organize  the  significant  features  of  these 
data  into  the  proto-type  of  an  outer  belt  dynamic  model.  As  for  the  validity  of  this  model, 
its  tests  against  data  will  be  the  final  arbiter. 

Meanwhile,  our  preliminary  work  summarized  here  needs  to  be  completed  in  the  fol¬ 
lowing  principal  areas: 

a.  Completion  of  diffusion  work  throughout  the  data  of  the  event. 

b.  Examine  the  effects  of  outer  magnetosphere  scattering  upon  fluxes  and  spectra. 

c.  Incorporate  usage  of  more  realistic  geomagnetic  fields  that  include  tail  and  better- 
modeled  magnetopauses. 

d.  Compile  a  storm  data  bcise,  using  the  August  1990  storm  as  benchmark. 

e.  Determine  short-term  storm  response  of  the  outer  belts  from  by  phenomenological 
and  theoreticcd  points  of  view. 

f.  Initiate  final  phase  of  outer  belt  dynamic  modeling  as  specified  in  the  proposal. 


11 


7.  Publications 


1.  Y.  T.  Chiu,  M.  A.  Rinaldi,  and  R.  W.  Nightingale,  Toward  dynamic  modeling  of 
the  outer  electron  radiation  belt,  J.  Geophys.  Res.,  95,  12069,  1990. 

2.  Y.  T.  Chiu,  W.  E.  Francis,  and  M.  Schulz,  An  exact  potential  model  of  the  outer 
magnetosphere,  J.  Geophys.  Res.,  in  preparation,  1991. 

3.  Y.  T.  Chiu,  M.  A.  Rinaldi,  W.  E.  Francis,  R.  M.  Robinson,  and  R.  W.  Nightingale. 
Outer  magnetosphere  and  radiation  belt  particle  dynamics  in  the  storm  of  August  26, 
1990,  Invited  paper,  American  Geophysical  Union  Spring  Meeting,  June  1,  1991. 

4.  M.  Schulz,  Y.  T.  Chiu,  and  W.  E.  Francis,  A  model  of  the  outer  magnetosphere. 
Invited  paper,  lUGG  General  Assembly,  August  12,  Vienna,  Austria. 


12 


MAGNETOSPHERE  IN  PARABOLOIDAL  HARMONICS 


o  O 

II  1> 

O  I'l 

>  ip:i 


o 

I — 
< 
ZD 

o 

LJJ 

LU 

O 

< 

_J 

Q. 

< 


c 


c 

c 

+ 

c 


c 

03 


di 

c 


.4" 

II 

HI 

O) 

3 

< 

CL 

o 

H 

Li] 

Z 

CD 

< 


O 


o 

I — 


o 

I — 
ZD 

_J 

o 

CO 


O 
C/D  O 


S  N 

II 

LJJ 

CD  O 


c  ^ 


z 

O 

O 

LU 

CC 

DC 

O 

DC 

LU 

I- 

Z 


o 


00 


o 

O 

ZL 

O  >  O 

%  Q  % 
B  cn  B 

O 

I- 


ZL 
h- 
< 
LU 
Z 

J  O 

03  h- 
LU 
Z 
CD 
< 


O 

ii 

=1 


z  •' 
O 

CD  O 
LU 


CC 

CC 

o 

CC 

LU 

1- 

X 

LU 


II 

=1 


Figure  1 
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BOW  SHOCK  &  OPEN  WAGNETOPAUSE  FIELD  LINE  TRACE 
TE=  8-APR-91.  SOLAR  WIND  (V  =1403  Km/s,  N  =  5  l/cm3) 

EXTERIOR  (Bx.By.Bz)  IS  (-10.0,  0.0,  0.0)  n.T,  E  =  0.20 

MAGNETOPAUSE  VALUES  NT.XM.ZM  &  A  ARE:  20,  6.3,10.0,  10. 
FIELD  LINE  IN  (X,Z)  PLANE,  DIPOLE  TILT=  0.0  DEG. 


NORTH 


CLOSED  "SSM"  FIELD  LINE  TRACES  NEAR  MAGNEiOPAUSE 
DATE=14-JUN-91,  solar  WIND  (V  =  254  Km/s,  N  =  10  Vcm3) 
EXTERIOR  (Dx.By.Br)  IS  (  0.0,  0.0,  10.0)  nT,  E  =  0.00 

MAGNETOPAUSE  VALUES  NT.XS.XM  &  XLIM  ARE;  15.12.6,10.0,-70 
FIELD  LINE  IN  (X,Z)  PLANE,  DIPOLE  TILT=  0.0  DEG. 
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Orbit  «  77o  E  -  800  koV  n  -  0  _ Orbit  =  77o  E  -  800  keV 


400  l<eV  n  =  0  _ Orbit  =  82o  E  =  400  keV 
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Figure  10 
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Spectrum  Before  Dropout 


Energy  (keV) 


.  Energy  Spectrum  During  Dropout 


00'99:W9l  Odd'AOy3N3Xmj  c: 


P3 


Figure  15 
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Spectrum  After  Dropout 


Energy  (keV) 


NORTH 


BEFORE 

DATE=26-AUG-90,  SOLAR  WIND  (V  =  701  Km/s,  N  =  5  /cm3) 

EXTERIOR  (Bx.By.Bz)  IS  (-10.9,  11.9,-11.8)  nT,  E  =  0.2C 
MAGNETOPAUSE  VALUES  NT,XM,ZM  &  A  ARE:  20,  8.0,  9.6,  20. 


20  10  0  -10  -20  -30 


SUN  <--  X 

Figure  17 
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NORTH 


ARER 


DATE=26-AUG-90.  SOLAR  WIND  (V  =1663  Km/s,  N  =  5  1/cm3) 

EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT,  E  =  0.20 

MAGNETOPAUSE  VALUES  NT.XM.ZM  <Sc  A  ARE:  20,  6.0,  7.2,  20. 


SUN  < —  X 


Figure  18 
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DUSK< - Y - >DAWN 


BOW  SHOCK  &  OPEN  M^GNETOPAUSE  FIELD  LINE  TRACE 
DATE  =  26-FEB-91.  SOLAR  WIND  (V  =1663  Km/s,  N  =  5  1/cm3) 

EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.S)  nT,  E  =  0.20 

MAGNETOPAUSE  VALUES  NT,XM,ZM  &  A  ARE:  20,  6.0,  7.2,  20. 
FIELD  LINE  IN  (X,Y)  PUNE,  DIPOLE  TILT=  0.0  DEG. 


SUN  <--  X 


Figure  19 
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RTH 


BOW  SHOCK  &  OPEN  MAGNETOPAUSE  FIELD  LINE  TRACE 
DATE=26-FEB-91.  SOLAR  WIND  (V  =1663  Km/s,  N  =  5  1/cm3) 

EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT,  E  =  0.20 

MAGNETOPAUSE  VALUES  NT,XM.ZM  &  A  ARE:  20.  6.0,  7.2,  20. 


FIELD  LINE  IN  (Y,Z)  PLANE,  DIPOLE  TILT=  0.0  DEG. 


DAWN< — Y — >DUSK 


Figure  20 
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DUSK< - Y - >DAWN 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.6,  0.0,  -4.0),  P.A.  =  90.2  deg 

DATE=14-FEB-91.  solar  WIND  (V=  701  km/s,  N=  5  1 /cm3) 
EXTERIOR  (8x,By,Bz)  IS  (-10.9,  11.9,-11.8)  nT,  E=  0.20 
Bext  =  20.02"[l.+0.00*SIN(  1.0»PHI+  0.)  +  0.00*SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM  =  2M  ARE:  8.0,  9.6 

ELECTRON.  Q  =  -1.  ENERGY=  1 0.004  Mev 


SUN  < —  X 


Figure  21 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.6,  0.0,  -4.0),  P.A.  =  90.2  deg 

DATE=14-FEB-91.  solar  .VINO  (V=  701  km/s,  N=  5  1/cti3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9.  11.9,-11.8)  nT.  E=  0.20 
Bext=20.02»[l.+0.00*SIN(  1.0»PHI+  0.)+0.00«SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  8.0,  9.6 

ELECTRON.  Q=-1.  ENERGY=  10.004  Mev 
TRAJECTORY  ON  (Y.Z)  PLANE.  DIPOLE  TlLT=  0.0  DEG. 


DAWN<--Y — >DUSK 


Figure  23 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.6,  0.0.  -4.0).  P.A.  =  90.1  deg 

DATE=14-FEB-91,  SOLAR  WIND  (V=1663  km/s.  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6.  3.1.  3.8)  nT,  E=  0.20 

Pext=  5.00*[1.+0.00*SIN(  1.0»PHI+  0.)+0.00*SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6  0.  7.2 

ELECTRON.  Q=-1.  ENERGY=  10.004  Mev 
TRAJECTORY  ON  (X.Y)  PLANE.  DIPOLE  TILT=  0.0  DEG. 


0 


NORTH 


(DEPT.  91-20),  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.6,  0.0,  -4.0),  P.A.  =  90.1  deg 

DAT£=14-FEB-91,  sour  WIND  (V=1663  km/s.  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT,  E=  0.20 

BexU  5.00»[1.+0.00»SIN(  1.0»''H1+  0.)  +  0.00*SIN(  1 .0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  Sc  YM=2M  ARE:  6.0,  7.2 

ELECTRON,  Q=-1,  ENERGY=  10.004  Mev 


TRAJECTORY  ON  (X.Z)  PUNE,  DIPOLE  TILT=  0.0  DEG. 


SUN  < —  X 


Figure  25 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.6,  0.0.  -4.0),  P.A.  =  90.1  deg 

DATE=14-FEB-91.  solar  WIND  (V=1663  km/s.  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6.  3.1.  3.8)  nT,  E=  0.20 

Bext=  5.00*[1.+0.00*SIN(  1.0*PHI+  0.)+0.00*SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6.0.  7.2 

ELECTRON.  Q=-1.  ENERGY=  1 0.004  Mev 


TRAJECTORY  ON  (Y.Z)  PLANE.  DIPOLE  TILT=  0.0  DEG. 


DAWN< — Y-->DUSK 


Figure  26 
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DUSK< - Y - >DAWN 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.0,  -6.3,  0.0),  P.A.  =  90.8  deg 

DATE=15-FEB-91.  SOLAR  WIND  (V=  701  km/s,  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9.  11.9.-11.8)  nT.  E=  0.20 
Bext=20.02*[l.+0.00*SIN(  1 .0*PHI+  0.)+0.00*SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &.  YM  =  ZM  ARE:  8.0.  9.6 

ELECTRON.  Q=-1.  ENERGY=  1 0.004  Mev 


TRAJECTORY  ON  (X.Y)  PLANE.  DIPOLE  TILT=  0.0  DEG. 


SUN  < - X 


Figure  27 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.0,  -6.3,  0.0).  P.A.  =  90.8  deg 

DATE  =  15-FEB-91.  SOLAR  WIND  (V=  701  km/s,  N=  5  l/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9.  11.9.-11.8)  nT,  E=  0.20 
Bext=20.02*[l.+0.00»SIN(  1.0*PHI+  0.)+0.00*SlN(  1.C*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  8.0.  9.6 

ELECTRON.  Q=-1.  ENERGY=  10.004  Mev 


TRAJECTORY  ON  (X.Z)  PLANE.  DIPOLE  TlLT=  0.0  DEG. 


SUN  < - X 


Figure  28 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.0,  -6.3,  0.0),  P.A.  =  90.8  deg 

DATE=15-FEB-91.  SOLAR  WIND  (V=  701  km/s,  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9,  11.9,-11.8)  nT.  E=  0.20 
Bext  =  20.02*[l.+0.00*SIN(  1.0*PHI+  0.)+0.00»SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM  =  ZM  ARE;  8.0,  9.6 

ELECTRON.  Q=-1,  ENERGY=  10.004  Mev 


TRAJECTORY  ON  (Y.Z)  PLANE.  DIPOLE  TILT=  0.0  DEG. 


DAWN<--Y-->DUSK 


Figure  29 
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DUSK< - Y - >DAWN 


AFTER 

STARTED  AT  (  1,0.  -6.3,  0.0),  P.A.  =  90.1  deg 

DATE  =  26-AUG-90.  SOLAR  WIND  (V=  1 663  km/s,  N=  5  1 /cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT,  E=  0.20 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM  =  2M  ARE:  6.0,  7.2 

ELECTRON,  0=-1,  ENERGY=  10.004  Mev 


SUN  < —  X 


Figure  30 
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NORTH 


(DEPT.  91-20),  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.0,  -6.3.  0.0),  P.A.  =  90.1  deg 

DATE  =  15-FEB-91.  SOLAR  WIND  (V=1663  km/s,  N=  5  "l /cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT,  E=  0.20 

Bext=  5.00*[1.  +  0.00»SIN(  1.0»PHI+  0.)+0.00»SIN(  1 .0»PHH-  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6.0,  7.2 

ELECTRON.  Q=-1,  ENERGY=  1 0.004  Mev 


TRAJECTORY  ON  (X.Z)  PLANE.  DIPOLE  TILT=  0.0  DEG. 


SUN  < —  X 


Figure  31 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  1.0,  -6.3,  0.0),  P.A.  =  90.1  deg 

DATE=15-FEB-91.  SOLAR  WIND  (V=1663  km/s.  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT,  E=  0.20 

Bext=  5.00*[1.+0.00»SIN(  1.0*PHI+  0.)+0.00*SIN(  1 .0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6.0,  7.2 

ELECTRON.  Q=-1.  ENERGY=  10.004  Mev 


TRAJECTORY  ON  (Y.2)  PLANE,  DIPOLE  TILT=  0.0  DEG. 


DAWN< — Y-->DUSK 


Figure  32 
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DUSK< - Y - >DAWN 


BEFORE 

STARTED  AT  (  7.0,  3.0,  0.0).  P.A.  =  90.9  deg 

DATE  =  26-AUG-90.  SOLAR  WIND  (V=  701  km/s,  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9,  11.9.-11.8)  nT,  E=  0.20 
MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM  =  ZM  ARE:  8.0,  9.6 

ION  Q=  1.  MASS=  1.  AMU.  ENERGY=  1.00  Mev 
TRAJECTORY  ON  (X.Y)  PLANE,  DIPOLE  TILT=  0.0  DEG. 
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Figure  33 
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(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  7.0,  3.0,  0.0),  P.A.  =  90.9  deg 

DATE=  8-FEB-91,  SOLAR  WIND  (V=  701  km/s.  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9,  11.9,-11.8)  nT,  E=  0.20 
Bext=20.02»[l.+0.00*SIN(  1.0»PHI+  0.)+0.00»SIN(  1.0»PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  YM=ZM  ARE:  8.0,  9.6 


ION  0=  1,  MASS=  1.  AMU,  ENERGY=  1.00  Mev 
TRAJECTORY  ON  (X.Z)  PLANE.  DIPOLE  TILT=  0.0  DEG. 


SUN  < —  X 


Figure  34 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  7.0,  3.0,  0.0),  P.A.  =  90.9  deg 

DATE=  8-FEB-91.  SOLAR  WIND  (V=  701  km/s,  N=  5  1 /cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (-10.9,  11.9,-11.8)  nT,  E=  0.20 
Bext  =  20.02*[l.+0.00*SIN(  1.0*PHI+  0.)+0.00*SIN(  1.0*PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM  =  2M  ARE:  8.0,  9.6 

ION  0=  1.  MASS^  1.  AMU.  ENERGY=  1.00  Mev 


DAWN<--Y-->DUSK 


Figure  3B 
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DUSK< - Y - >DAWN 


AFTER 

STARTED  AT  (  -4.0,  5.9,  -0.1),  P.A.  =117.2  deg 

DATE=26-AUG-90,  SOLAR  WIND  (V=1663  km/s,  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6,  3.1,  3.8)  nT.  E=  0.20 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6.0,  7.2 


ION  Q=  1,  MASS=  1.  AMU,  ENERGY=  1.00  Mev 


SUN  <--  X 


Figure  36 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  -4.0,  5.9.  -0.1).  P.A.  =117.2  deg 

DATE=1 1 -FEB-91.  SOLAR  WIND  (V=1663  km/s.  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6.  3.1,  3.8)  nT.  E=  0.20 

Bext=  5.00»[1.+0.00»SIN(  1.0-PHI+  0.)+0.00*SIN(  1.0»PHI+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6.0.  7.2 


ION  0=  1.  MASS=  1.  AMU.  ENERGY=  1.00  Mev 
trajectory  on  (X.Z)  plane.  DIPOLE  TlLT=  0.0  DEG. 


SUN  <--  X 


Figure  37 
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NORTH 


(DEPT.  91-20).  LPARL  PARTICLE  TRAJECTORY  CODE 
STARTED  AT  (  -4.0,  5.9,  -0.1).  P.A.  =117.2  deg 

DATE=1  1 -FEB-91,  SOLAR  WIND  (V=1663  km/s,  N=  5  1/cm3) 
EXTERIOR  (Bx.By.Bz)  IS  (  0.6.  3.1,  3.8)  nT.  E=  0.20 

Bext=  5.00»[1.+0.00*SIN(  1.0*PHI+  O.)-i-C.O0*SIN(  1.0-*''’'l+  0.)] 

MAGNETOPAUSE  BOUNDRY  VALUES  XM  &  YM=ZM  ARE:  6.0,  7.2 


-50  5  10 

DAWN<  —  Y-->DU5K 


Figure  38 
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Pilch-Angle  (deg) 


